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The relationship between the magnetic field and the circular polarization of astro- 
physical maser radiation due to the Zeeman effect under idealized conditions is inves- 
tigated when the Zeeman splitting is much smaller than the spectral linebreadth and 
when radiative saturation is significant. The description of the circular polarization, as 
well as inferences about the magnetic field from the observations, are clearest when the 
rate for stimulated emission is much less than the Zeeman splitting. The calculations 
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ON . here are performed in this regime, which is relevant for some (if not most) observations 

of astrophysical masers. We demonstrate that Stokes- y is proportional to the Zeeman 
splitting and the fractional linear polarization is independent of the Zeeman splitting 
when the ratio of the Zeeman splitting to the spectral line breadth is small — less than 
about 0.1. In contrast to its behavior for ordinary spectral lines, the circular polarization 
Q H i for masers that are at least partially saturated does not decrease with increasing angle 

between the magnetic field and the line-of-sight until they are nearly perpendicular. 
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The circular polarization of maser radiation where the Zeeman splitting is much smaller than 
the spectral linebreadth is utilized in efforts to obtain information about the magnetic field in 
various astronomical contexts (e.g., for recent work, Kemball & Diamond 1997; Sarma, Troland, 
& Romney 2001; Vlemmings, Diamond, & van Langevelde 2001; Yusef-Zadeh et al. 1999). For a 
two-level, non-masing ( "thermal" ) spectral line, a simple relationship exists 

V/(dI/dv) =pB cos 6 (1) 

that involves Stokes-F, the derivative of the intensity I with respect to Doppler velocity v, the 
magnetic field B, and the angle 9 between the magnetic field and the line-of-sight (p is a constant 
of the molecular physics for the specific transition). When maser saturation is unimportant, the 
relationship can readily be seen as applicable for masers as well (e.g., Fiebig & Gusten 1989). 



Though maser polarization in the presence of saturation was the focus of the classic investigation 
of Goldreich, Keeley & Kwan (1973; hereafter GKK), they provide no guidance about the circular 
polarization for weak Zeeman splitting. GKK only consider line center, where the circular polar- 
ization is zero in the weak splitting regime. For a summary of the theory of maser polarization, see 
Watson (2001). 

The GKK idealization is a uniformly pumped, linear maser in steady-state. The seed radiation 
is external, weak continuum radiation and the masing involves an angular momentum J = 1 — 
transition in the presence of a constant magnetic field. Even though the actual circumstances 
may be more complicated, the idealization and solutions of GKK have served as the basis for 
discussions of maser polarization in general. In this Letter, we partially remedy the omission of 
GKK for weak splitting by calculating Stokes- F for this basic idealization in the limit that the rate 
R for stimulated emission is much less than the Zeeman splitting gVt (where 17 = eB/m e c and g is 
the Lande g- factor). This limit certainly is applicable when gQ is a noticeable fraction (greater than 
a few percent) of the linewidth as occurs for the 1720 MHz OH masers. It probably also is relevant 
for the 22 GHz water masers for which g£l ~ 10 4 B{G)s~ 1 for the F = 7 — 6 transition and provides 
at least a basis for discussion for the SiO masers where gQ, ~ 10 3 -B(G)s -1 . As a benchmark, the 
measure of the degree of saturation is the ratio of R to the loss rate T which ordinarily is taken 
to be about 1 s _1 and 5 s _1 for the H2O and SiO masers, respectively, and somewhat less for 
the 1720 MHz masers. We limit our attention to g£l 3> R because here the modifications to the 
standard Zeeman effect are moderate and the magnetic field can still be inferred from Stokes- V with 
some confidence. Both the conceptualization and calculation of Stokes- F are considerably simpler 
in this regime because only "ordinary rate equations" and "ordinary molecular populations" of 
the magnetic substates are involved. In contrast, when gQ ^> R is not satisfied, Stokes-V^ can 
be much larger (sometimes referred to as "non-Zeeman" effects) than expected from the standard 
Zeeman relationship (Nedoluha &; Watson 1994). It is unclear whether useful information about the 
magnetic field can be inferred from the observations in this regime. The more involved methods 
of the quantum mechanical density matrix must be utilized in this regime to incorporate the 
correlations between the magnetic substates. Large, non-Zeeman circular polarization can also be 
generated by effects that would not be present in the GKK idealization for masing, such as changes 
in the direction of the magnetic field within the masing region (Wiebe & Watson 1998). 

Calculations of the type being presented here have previously been performed (Nedoluha & 
Watson 1992), but only for a quite limited range of saturation and angles 6, and only for the 
molecular parameters that are specific to the 22 GHz masing transition. That is, for the high 
angular momenta (F = 7, 6, 5, & 4) and when the values of the magnetic moments are different for 
the upper and lower molecular states. These calculations did demonstrate that maser saturation 
alters equation (1) and that its effect is to increase Stokes-y much more at large angles than 
at small angles relative to that expected from equation (1). The latter is especially significant 
because lines-of-sight to masers are often likely to be nearly perpendicular to the magnetic fields 
(e.g., toward circumstellar masing rings, toward the centers of accretion disks viewed edge-on). 



The intent of this Letter is to provide a comprehensive description of the relationship between the 
circular polarization and the magnetic field within the limitations outlined above, and to do this 
for the idealized masing conditions treated by GKK which have served as a basis for discussing 
maser polarization. 



2. CALCULATIONS 

The rate equations for the normalized (dimensionless) population differences n+,ri-, and uq 
between the magnetic substates m = +1, -1, and of the upper (angular momentum J = 1) energy 
level and the m = substate of the lower (J = 0) level, respectively, are functions of the molecular 
velocity v and can be expressed by (e.g., Wallin & Watson 1995; also GKK) 

1 = {l + 2R + )n + + R n + R-n_ (2) 

and the two related equations that are obtained by exchanging the ± and indexes in equation 
(2). "Phenomenological" pumping and loss rates "A" and 'T", which are standard in maser theory, 
have been utilized and incorporated into the normalizations of the populations to obtain equation 
(2). For a unidirectional linear maser, the normalized rates for stimulated emission in equation (2) 
are 

R± = I± (l + cos 2 6) + Q±sin 2 6 ± 2V±cos<9 (3) 

and 

R = 2(I -Q )sm 2 e (4) 

where the subscripts ± and indicate that the intensities are to be evaluated at the frequencies 
(angular) u± and ujq for which 

{l-v/c)u± = u R ±gn/2 (5) 

and 

(1 - v/c)uj = w fi (6) 

Here, lor is the resonance frequency of the transition and ±g£l/2 are the Zeeman shifts of the m = 
±1 magnetic substates. The dimensionless Stokes intensities (I, Q, V) are actual intensities divided 
by the characteristic saturation intensity Is = (8?zo; 3 r /3irc 2 Ae) where Ae is the Einstein A-value 



for the transition. Note that whereas R± and i?o are dimensionless, R itself is not normalized and 
retains its usual dimensions (time -1 ) in our discussions. The radiative transfer equations 

dl/ds =AI + BQ + CV (7) 

dQ/ds = AQ + BI (8) 

dV/ds = AV + CI (9) 

are at a specific frequency uj and can be solved numerically as a function of the normalized distance 
s with the expressions 

A= (1 + cos 2 6)(f + n + + f_n^) + 2f n sm 2 6 (10) 

B = sin 2 9{f + n + + /_n_ - 2/ n ) (11) 

C = 2cos0{f+n + - /_n_) (12) 

where the /'s are Maxwellian distributions for the component of the velocity that is along the path 
of the radiation. In these equations, f±,n±,fo and no are evaluated at velocities v± and vq given 
by (1 — v±/c)uj = ujr ± gQ/2 and (1 — vq/c)uj = ujr. 

Starting from the initial conditions /< 1 and Q = = V, the radiative transfer equations 
are integrated numerically at a sufficient number of frequencies u to delineate the profile of the 
spectral line. We emphasize that the foregoing equations are valid only when gQ S> R- They are 
the same as would describe the polarization of ordinary thermal radiation — the only differences 
being the sign of the net pumping rate which is reflected in the sign of the "1" on the left-hand-side 
of equation (2) and the ignoring of spontaneous emission. 

The main results are presented in the Figure 1 in terms of the ratios V/(pBdI/dv) and Q/I 
computed for a value of gQ that is much less than the spectral linebreadth and at a frequency 
that is essentially (though not exactly) at line center. These ratios are useful because, in the 
weak splittng regime and consistent with equation (1), they tend to be independent of g£l and 
constant with frequency across the spectral line. At frequencies where the bulk of the intensity 
occurs, we have confirmed that the ratios in Figure 1 are independent of frequency to a good 
approximation in the weak splitting regime as delineated in the following paragraph. Exceptions 
to this generalization do occur at the angles for which the deviations from equation (1) are greatest 



and also the intensities are in the neighborhood of 1=1. Where there are variations in V/(pBdI/dv) 
and Q/I with frequency, these variations are modest and are smallest near line center. Note that 
the properties of the radiation at different frequencies within the spectral line are not completely 
independent because of the coupling that is indicated in equations (5) and (6), and in the analogous 
expressions for v± and vq. In detail, the rate for stimulated emission depends on the polarization 
of the radiation, on the angle 9, and is not exactly the same for all of the magnetic substates. It 
nevertheless is useful to utilize / (which is normalized by the characteristic saturation intensity) as 
the measure of the degree of saturation. That is, / ~ R/T. 

The goal here is limited to providing information on the polarization characteristics for the limit 
of weak Zeeman splitting — that is, for small enough gQ that the ratios in Figure 1 are independent 
of gQ to a good approximation for the bulk of the radiation within the spectral line. The precise 
accuracy of this approximation for a specific gU depends upon angle 6, the degree of saturation, the 
Doppler velocity within the spectral line, and is not exactly the same for the linear and the circular 
polarization. Nevertheless, some useful generalizations can be made. We compute polarizations for 
gn as large as one-fifth of the FWHM thermal Doppler breadth Au t (Au t = 2Aoj R [kT/Mc 2 ] l/2 ). 
As long as the degree of saturation is greater than one, the ratios presented in the Figure 1 are 
independent of gVt to an accuracy of a few percent. Likewise, as long as gQ/AuJt < 0.05, these 
ratios also are independent of gQ, to a similar accuracy regardless of the degree of saturation. For 
lower saturations — I down to 0.01 — and 0.2 > gtt/Aut > 0.05 , the ratio V/(pBdI/dv) is still 
independent of g$l to within about twenty percent accuracy. The main deviations occur at angles 
less than about forty-five degrees where the effects of saturation are the least. 

The main effect of saturation on the circular polarization is readily evident in the Figure 1. 
Instead of decreasing as cos 6 when the angle 6 increases from zero, Stokes- F tends to increase 
until a relatively large angle is reached whose value is determined by the exact degree of saturation. 
At 6=0, Stokes-V becomes equal to its value for unsaturated masing. The nearly discontinuous 
change reflects the change in the molecular populations that is associated with the similarly rapid 
variation of the linear polarization as 6 approaches zero in the GKK theory. The variation of the 
linear polarization with saturation can be sensitive to the angular momentum of the molecular 
states (e.g., Deguchi &: Watson 1990). We have thus performed similar calculations for a J = 2 — 1 
masing transition, as well, with the idealization that g is the same for the upper and lower energy 
levels (hence, p = ge/2m e UR is the same for both transitions). These also are presented in Figure 
1, along with additional results from a further calculation in which the J = 1 — masing is treated 
as bidirectional. The optical depths are the same for radiation propagating in opposite directions 
in a linear maser. If the external seed radiation is similar at both ends of the maser (or if the 
seed radiation is due to spontaneous emission), the radiation propagating in the two directions will 
also be similar and the maser will be "bidirectional" . It is known that the variation of the linear 
polarization with saturation can be different for uni- and bidirectional masers (Western &: Watson 
1984). However, the comparisons in Figure 1 demonstrate that the differences for the circular 
polarization are small for both of these modifications which we have considered beyond the basic 



unidirectional J = 1 — maser . 

The fractional linear polarization in Figure 1 is identical with that obtained in previous cal- 
culations. We present it here in the same format as the circular polarization for convenience in 
relating the two polarizations. At high saturation, the fractional linear polarizations in our calcula- 
tions (see also Western & Watson 1984) do reach those of GKK — which are obtained by GKK only 
for the limit of high saturation. The fractional linear polarization is seen in Figure 1 to increase 
more rapidly as a function of saturation for J = 1 — than for J = 2 — 1 masers, and more rapidly 
for bidirectional than for unidirectional masers. As found previously, however, to reach the very 
highest Q/I given by GKK, the degree of saturation must be implausibly high for J = 2 — 1 masers. 
Anisotropic (or "m-dependent") pumping by starlight is thus the most likely cause for the highest 
fractional linear polarizations of the SiO and perhaps other masers (Western & Watson 1983). 

Unfortunately, the degree of radiative saturation of astrophysical masers is a longstanding 
uncertainty. In addition to the surface brightness of the maser, the saturation depends on the angle 
into which the radiation is beamed — a quantity for which direct estimates ordinarily are unreliable 
(e.g., Watson & Wyld 2000). Note that the "intensities" of the linear maser formulation (by us 
and by GKK) actually are "mean intensities" x4tt. Hence, intensities from the observations must 
be multiplied by the solid angle of the beam (as well as divided by 1$) in order to relate them to 
the / in the Figures. A comparison of the observed linear polarization with that in Figure 1 can be 
helpful to restrict the degree of saturation if there are no contaminating effects such as anisotropic 
pumping, multiple components, or Faraday rotation. The narrowing of the spectral line during 
unsaturated maser amplification, followed by rebroadening of the line when the maser becomes 
saturated, can also be used to obtain an indication of saturation (e.g., Nedoluha & Watson 1991) 
if there are no velocity gradients or multiple components within the maser. A consideration of the 
linebreadths of certain prominent 22 GHz water masers indicates that the influence of saturation 
on Stokes- V probably is modest for these masers (Nedoluha & Watson 1992). The linebreadth in 
our calculations is given in Figure 2 as a function of /. It does depend somewhat on the intensity of 
the incident continuum seed radiation. The computations in Figure 1 are performed for an incident 
intensity / = 10~ 9 which we believe to be representative. We have also performed computations 
analogous to those in Figure 1 when the incident intensity I is 10~ 5 as might occur when the masing 
gas is amplifying a strong continuum source. The resulting Q/I are essentially identical to those 
in Figure 1. The Stokes- U also are quite similar to those in Figure 1, except near the peaks of the 
curves for the the largest / (> 10 2 ) where they are smaller by 15-20%. 

In summary, when the observed masers are believed to be at least somewhat saturated, but 
there is no good information about the angle 9 nor about the exact degree of saturation, simply 
removing the cosO in equation (1) would seem to provide the best way at present to infer mag- 
netic field strengths from the observed Stokes-U in the weak splitting regime when "non-Zeeman 
effects" can be ignored. Saturation with / > 10 2 (and probably even I > 10) seems unlikely for 
astrophysical masers. In contrast to the linear polarization, the circular polarization is relatively 
insensitive to the angular momentum of the molecular states. We emphasize that our results are 



applicable in detail only for the idealized masing conditions on which the calculations are based 
(see Introduction). In addition to non-Zeeman effects, velocity gradients, anisotropic pumping, 
and multiple hyper fine components may be present, but are not considered here. For example, 
the 22 GHz masing transition of water probably consists of multiple hyperfine components so that 
equation (1) is unlikely to be directly applicable for it (but see Nedoluha &; Watson 1992). 
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Fig. 1. — Circular and linear polarization of maser radiation as a function of the cosine of the angle 
9 between the direction of the magnetic field and the line of sight. Polarizations are presented 
for J = 1 — unidirectional and bidirectional masers, and for J = 2 — 1 unidirectional masers, in 
separate panels as indicated by the label in each panel. The curves are labeled by the log 10 of the 
intensity I which is essentially the degree of saturation. Curves (some of which overlap) are plotted 
for I = 10~ 2 , 10 _1 , 1, 3, 10, 10 2 , 10 3 , and 10 4 . 

(upper three panels) The circular polarization as measured by the magnitude of V/(pBdI/dv) 
which is equal to cos# in the unsaturated limit. 

(lower three panels) The fractional linear polarization. Stokes-Q is positive when the linear po- 
larization is perpendicular to the direction of the magnetic field projected onto the plane of the 
sky. 

Fig. 2. — The ratio of the spectral linebreadth (FWHM) of the maser radiation to the Doppler 
breadth Atot (FWHM) of the thermal velocities of the masing molecules, as a function of the log 10 
of the intensity I. Linebreadths are shown for J = 1 — unidirectional (solid lines), J = 1 — 
bidirectional (dotted lines), J = 2 — 1 unidirectional (dashed lines) masers with intensities for the 
incident continuum radiation of / = 10~ 5 , as well as for I = 1CP 9 which we consider to be most 
representative. 
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